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Low temperature specific heat has been measured and extensively analyzed on a series of 
La2-xSrxCu04, single crystals from underdoped to overdoped regime. From these data the quasi- 
particle density of states (DOS) in the mixed state is derived and compared to the predicted scaling 
law Cvoi/TvTl = f{T/y/lI) of d-wave superconductivity. It is found that the scaling law can be 
nicely followed by the optimally doped sample (x=0.15) in quite wide region of (T/y/lI < 8K/Vt). 
However, the region for this scaling becomes smaller and smaller towards more underdoped region: 
a clear trend can be seen for samples from x= 0.15 to 0.069. Therefore, generally speaking, the 
scaling quality becomes worse on the underdoped samples in terms of scalable region of T/y/TI. 
This feature in the underdoped region is explained as due to the low energy excitations from a sec- 
ond order (for example, anti-ferromagnetic correlation, d-density wave, spin density wave or charge 
density wave order) that may co-exist or compete with superconductivity. Surprisingly, deviations 
from the d-wave scaling law have also been found for the overdoped sample (x=0.22). While the 
scaling law is reconciled for the overdoped sample when the core size effect is taken into account. An 
important discovery of present work is that the zero-temperature data follow the Volovik's relation 
A7(r = 0) = Ay/ll quite well for all samples investigated here although the applicability of the 
d-wave scaling law to the data at finite temperatures varies with doped hole concentration. Finally 
we present the doping dependence of some parameters, such as, the residual linear term 70, the a 
value, etc. It is suggested that the residual linear term {'yoT) of the electronic specific heat observed 
in all cuprate superconductors is probably due to the inhomogeneity, either chemical or electronic in 
origin. The field induced reduction of the specific heat in the mixed state is also reported. Finally 
implications on the electronic phase diagram are suggested. 

PACS numbers: 74.20.Rp, 74.25.Dw, 74.25.Fy, 74.72.Dn 



I. INTRODUCTION 

One of few points with consensus in the cuprate su- 
perconductors is the d.j;2_y2 pairing symmetry in hole 
doped region. This has been supported by tremen- 
dous experiments^ both from surface detectionSiiii^ 
and bulk measurementsiSiSiiSiiLiSiiiii. In a d-wave su- 
perconductor with hne nodes in the gap function, the 
quasiparticle density of states (DOS) N{E) rises linearly 
with energy at the Fermi level in zero field, i.e., N{E) oc 
\E — Ep\, resulting ini^ an electronic specific heat Ce = 
aT^, where a oc jn/Tc and 7„ is the specific heat which 
is proportional to the DOS at the Fermi level of the nor- 
mal state. In the mixed state with the field higher than 
a certain value, the DOS near the Fermi surface becomes 
finite, therefore the quadratic term Ce = ctT^ will be 
surpassed and substituted by both the localized excita- 
tions inside the vortex core and the de-localized excita- 
tions outside the core. Volovik^^ pointed out that for d- 
wave superconductors in the mixed state, supercurrents 
around a vortex core lead to a Doppler shift to the quasi- 
particle excitation spectrum, which affects strongly the 



low energy excitation around the nodes. It was shown 
that the contribution from the delocalized part (outside 
the core) will prevail over the core part and the specific 
heat is predicted to behave asi^^ Cyoi = k^fnT^/W/H^ 
with k in the order of unity. This prediction has been ver- 
ified by many measurements which were taken as the evi- 



dence for d-wave symmetrji 
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In the finite 



temperature and field region a scaling law is proposediS 
as 



aoi/TVH = fiT/VH) 



(1) 



with T/Vh < Tc/^jHc2{0). This scahng law can be 
further converted into the form of C-aoijU = 9{T /\/H) 
or Cyoi/T'^ = y{T/VH), here f{x), or g{x) or y{x) are 
unknown scaling functions. This scaling law has been 
proved in YBCCM^iS. and in LSCOiiiSJ^ii^. It remains 
however unclear whether this scaling law is still valid 
in the very overdoped region since the vortex core size 
^ grows up. In the underdoped region, inelastic neu- 
tron scattering reveals that an anti-ferromagnetic order 



emerges when the superconductivity is suppressedi2i2S. 
It is thus also interesting to check whether the d-wave 
scahng law proposed by Simon and Lee is applicable in 
underdoped regime. In addition, the Simon-Lee scaling 
law is in agreement with the calculations as proposed 
by Volovik and Kopnini^ti^ in two extreme conditions of 
temperature. In the low temperature limit the scaling 
law Cvoi/T'^ = y{T /\/H) becomes the Volovik's relation 
Cvoi — AT\/H. When the temperature is increased, an- 
other relation C oc aT^ + bH is reached. The boundary 
between these two regions is T/^/H = Tc/\/Tid2 accord- 
ing to Volovik and Kopnin^i. These theoretical models 
can be quantitatively tested by experiments on samples 
with different doping concentrations. 

Another important but controversial issue is the vor- 
tex core state in the cuprate superconductors. By solv- 
ing the mean-field Bogoliubov-de Gennes (BdG) equa- 
tion, theoretically it is suggested that a zero-bias con- 
ductance peak (ZBCP) exists in the vortex cor o^^i^? . 
However this is in sharp contrast with the experimental 
observationsS^iSSiSSiSLSSiS, mainly on optimally doped 
samples. The absence of a ZBCP within the vortex cores 
was attributed to the presence oiidxy or is components'^, 
or the competing orders (see later) . In this paper we show 
that the DOS due to vortex quasiparticle excitations de- 
viates from Simon-Lee scaling law for the overdoped sam- 
ple, but follows rather well with the optimally doped one. 
The deviations for the overdoped sample are found to be 
induced by the vortex core size effect. In the extremely 
underdoped region, it is found that the Simon-Lee scal- 
ing law fails, except for in very low temperature region. 
This can be understood as due to the competing order 
emerging within or nearby the vortex cores. 



II. EXPERIMENT 

The single crystals measured in this work were pre- 
pared by the travelling solvent floating-zone technique. 
Samples with seven different doping concentrations 
(p=0.063, 0.069, 0.075, 0.09, 0.11, 0.15, 0.22) have been 
investigated. The sample with p=0.15 and 0.22 are from 
GRIEF, and others are from NLSG(IOP). Part of the 
data for all samples will be presented, for example the 
field induced change of 7 at zero K, the residual lin- 
ear term 70 and the a value in the pure d-wave ex- 
pression Cdos — OiT'^ when H — Q (see later). How- 
ever for clarity we mainly show data and the analysis 
on three typical samples with x=0.22 (Tj. = 27. 4K, over- 
doped), x=0.15 {Tc = 36. IK, close to optimal doping 
point) and p=0.069 {T^ w 12K, underdoped, x=0.063 
originally) as characterized by AG susceptibility and DG 
magnetization (shown by the insets in Fig. 3, Fig. 8 and 
main panel of Fig. 12). The quality of our samples has 
also been characterized by x-ray diffraction patterns, and 
R{T) data showing a narrow transition ATc < 2 K. 
For some samples, the full width at the half maximum 
(FWHM) of the rocking curve of the (008) peak is only 



0.10°. The overdoped sample has a mass about 28.55 
mg and 3.66 x 2.3 x 0.577im'^ in dimension. The opti- 
mally doped sample weighs about 23.6 mg with dimen- 



sions of 3.1 X 3 X 0.5f 



For the underdoped sample 



with nominal concentration x=:0.063, before annealing, it 
has a superconducting transition temperature of about 
12 K and a mass of about 32.89mg and 3.75 x 2.75 x 
0.5?nTO'^ in dimensions. By fitting to the empirical rela- 
tion T^/T^^'''^ = 1 - 82.6(p - 0.16)2 with T™*^^ = 38X 
the maximum T,. at the optimal doping point p — 0.16, 
we estimate that the hole concentration of this sample is 
around p — 0.069. After annealing in flowing Ar gas for 
48h, the Tc drops down from about 12 K to 9 K indi- 
cating that the sample becomes more underdoped. Note 
that Tc=^ K is expected exactly by the empirical relation 
aX X = p — 0.063. 

The heat capacity presented here were taken with the 
relaxation method based on an Oxford cryogenic sys- 
tem Maglab. The sample is put onto a micro-chip on 
which there is a tiny Gernox temperature sensor, a film 
heater. The micro-chip together with the sample are 
hung up by golden wires in vacuum. These golden wires 
are the only thermal links between the micro-chip and 
the thermal sink whose temperature is well controlled. 
The temperature of the micro-chip is controlled by the 
on-board small film heater and measured by the on-board 
thermometer. When the temperature of the micro-chip 
is stable, a heating power with fixed current is sent to 
the film heater on the chip and the time dependence of 
the chip temperature is measured simultaneously. The 
change of temperature is fitted to an exponential rela- 
tion AT = ATo X [1 — ea::p(— i/r)], and heat capacity 
is determined by r = {C + Cadd)/^^, here C and Cadd 
are the heat capacity of the sample and addenda (includ- 
ing a small sapphire substrate, small printed film heater, 
tiny Gernox temperature sensor, 025 fj,m gold wire leads, 
Wakefield thermal conducting grease (about lOOfig)) re- 
spectively, where k^ is the thermal conductance between 
the chip and the thermal link. The value C'add has been 
measured and subtracted from the total heat capacity, 
thus C value reported here is only that from the sample. 
In Fig. 1 we present the temperature and field dependence 
of the heat capacity from the addenda and three typical 
samples. It is clear that the heat capacity of the addenda 
is much smaller than the value of the samples. In addi- 
tion, the data of C/T for the addenda extrapolates to 
zero at r = OK showing the only existence of phonon 
part. One can also see that almost no field dependence 
can be observed for the addenda. However, for all sam- 
ples, there is a clear finite intercept at T = OK which 
gives rise to a residual linear term 70. Meanwhile the 
field induced change can be easily observed for all sam- 
ples, even for the very underdoped sample. The inter- 
crossing of the data at H ^ and H = 12T at about 6K 
for the undedoped sample is understandable and will be 
discussed later. 

In all measurements done for the present work, the 
magnetic field H is always parallel to c-axis of single crys- 
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FIG. 1: Temperature dependence of the heat capacity from 
the addenda with 110 [ig Wakefield grease (bottom sohd line 
for H=0 and dashed line for H=12 T), and three typical sam- 
ples (filled symbols for H=0 and open symbols for H=12 T, 
lines are guides to the eye). 



tals, and the data are collected in the warming up pro- 
cess after it is cooled under a field (Field-Cooling pro- 
cess, abbreviated as FC hereafter). In the data treat- 
ment we use A7 = [C'/f||c — C//=o]/T instead of using 
A7 = ^H\\c — Chlc]/T. The latter may inevitably 
involve the unknown DOS contributions from another 
kind of vortices (for example, Josephson vortices) when 
iJ _L C. The field dependence of the Cernox thermometer 
has been calibrated well by Oxford before the shipment. 
The true temperature has been derived automatically by 
the software with a calibration table with magnetic fields 
at OT, IT, 2T, 4T, 8T and 12 T. The values at other 
fields are obtained also automatically by software by do- 
ing linear interpolation between two nearby fields. There- 
fore the readout from the machine gives directly the true 
temperature value with the field effect corrected. 



III. RESULTS AND DATA ANALYSIS 

A. Fitting to the zero field data 

Before showing the field induced change of the heat ca- 
pacity we present in Fig. 2 the temperature dependence of 
C /T for some samples at zero field. As mentioned previ- 
ously, for a d-wave superconductor in the superconduct- 
ing state, it is known that Coos = o/T"^ when H = 0. 
In addition, as observed in other cuprate superconduc- 
tors, the curve at zero field extrapolates to a finite value 
(70) at K instead of zero. This was interpreted as po- 
tential scattering effect due to small amount impurities 
or disorders2iffi. We will argue that this residual linear 
term may also reflect physics beyond the simple argu- 
ment of impurity scattering (see later) . As also observed 
by other groups for La — 214 system, the anomalous up- 
turn of C/T due to the Schottky anomaly of free spins 



FIG. 2: Temperature dependence of C/T for samples 
(p=0.069,0.075, 0.09, 0.11,0. 15 and 0.22) at zero field. The 
solid lines are fits to eq.2 and the parameters derived here 
are listed in Table-I. 
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is very weakiitiStiii^. This avoids the complexity in the 
data analysis. Together with the phonon contribution 
/3T^ + ST^, we have 



C{H = 0)/T ^jo + aT + (iT^ + 5T^ 



(2) 



Above equation is used to fit the data at if = for 
some samples. The fitting results are shown in Fig. 2 
and listed in Table-I, where the units for 70, a, (3 and 
5 are mJmol~^K~^, mJmol~^K~^, mJmol~^K~'^ and 
mJmol~^K~^ respectively. One can see that a decreases 
quickly towards underdoping, P (and thus the Debye tem- 
perature Qd) does not change too much with doping. 
The sudden drop of a at p = 0.11 may be induced by 
the well known 1/8 problem. The residual linear term 70 
increases rapidly towards underdoping, which will be dis- 
cussed later. The a values are also comparable to those 
found by other groupsi^*i^. 



B. Overdoped sample with x^O.22 

Fig. 3 shows C/T as a function of T^ at magnetic fields 
ranging from to 12 T for the overdoped sample. The 
separation between each field can be well determined. 




10 20 30 40 50 60 70 80 90 100 



T' {K') 










AAAAi ■ " V V V " O 

AAA V V 

■nn„nnnn°o° ^ A A a 

□ □ P Q Q Q D 



1 T 
4T 
8T 
12T 



< e 



o e Q 



10 



T(K) 



FIG. 3: Specific heat C/T vs. T^ of tlie overdoped sample 
(x — 0.22) at magnetic fields ranging from to 12 T. The 
inset shows the diamagnetic transition at around 27.4 K de- 
termined by the crossing point of the extrapolating line of the 
most steep part with the normal state background M = 0. 



In low temperature region the curves are rather linear 
showing that the major part is due to phonon contribu- 
tion Cph = (iT^ + 5T^ . It is known that the phonon 
part is independent on the magnetic field, this allows to 
remove the phonon contribution by subtracting the C/T 
at a certain field with that at zero field. The results af- 
ter the substraction are shown in Fig. 4. The subtracted 
values A7 = 7(H) - 7(0) = [C{T, H) - C{T, H = 0)]/T 
exhibit a rather linear T dependence in low temperature 
region. One can also see that the negative slope is actu- 
ally field dependent. In the following we will show that 
the field dependent slope of the linear part in low tem- 
perature region shown in Fig. 4 directly deviates from the 
Simon-Leei^ scaling law. 

According to Simon-Lee scaling law C^oilT\/li = 
f(T/\^), in low temperature region, the Volovik's rela- 
tion restores, thus one has Cyoi/T = A\fli and further 
has 



A7 = {C{H) - C(0)]/T = A^Jli - aT 



(3) 



This clearly shows that there is a negative slope for 
A7 vs. T, but the slope a is a constant. However, from 
Fig. 4 one can see that the slope changes slightly with the 
magnetic field H. This indicates that only the Volovik's 
relation is not enough to interpret the data. In the fol- 
lowing we will take both the core size effect and finite 
temperature effect into account. The former has not been 
considered in the original Simon-Lee scaling law since the 
size of the vortex core was thought to be small, and the 
contribution from that small region is negligible. If the 
vortex core size becomes bigger, this should be reconsid- 
ered when counting the DOS due to the Doppler shift 
effect. 

Let us first consider only the finite temperature effect. 
Suppose that we are in the crossover region between the 



FIG. 4: Temperature dependence of A7 = ^(H) — 7(0) = 
[C{H, T) - C(0, T)]/T of the overdoped sample (x = 0.22). A 
linear behavior is clearly seen in low temperature region with 
a field dependent slope, which is not in accord with the pro- 
posed scaling law by Simon and Lee (see text). The straight 
lines in low temperature region are guides to the eyes. From 
these lines one can determine the zero temperature intercept 
A7 and the slope dA'y/dT shown in Fig. 5 



low temperature limit and high temperature limit sug- 
gested by Volovik and Kopnin^''^, making Taylor's expan- 
sion to the right hand side of Simon-Lee scaling law leads 
to 



Cyoi = boH + biTVH + b2T^ + o{T^) 



(4) 



where 69 = because Cyoi /T should not diverge when 
T = and H ^ 0, 61 = A. Since o{T^) is very 
small in low temperature region, one thus has C^oi/T = 
bivH + 62T. Interestingly, one can see that this simple 
formula contains the results both in the low tempera- 
ture limit Cvoi = biTvHi^'i^'^ and high temperature 
limit^-'^-'^^ Cvoi = b2T^. This is not surprising since a 
scaling function should be more general and cover most 
possible cases. When Hd << H << Hc2, the total spe- 
cific heat contains four parts: Doppler shift term from 
the region outside the core Cvoh the inner vortex core 
term Ccore oc HT, the residual linear term 70T and the 
phonon term Cph- Here it is assumed that the heat ca- 
pacity contributed by the core region is equal for each 
vortex and independent on the external magnetic field, 
thus Ccore depends only on the vortex density which is 
proportional to H. The local DOS measured by STMS& 
revealed that the low energy DOS within the vortex core 
differs only slightly from the case for a d-wave supercon- 
ductivity (outside and far away from the vortex core). 
When changing the external magnetic field, the low en- 
ergy DOS within the vortex core is not expected to vary 
too much. At zero field, the total specific heat contains 
three parts: JqT and Cph, and a quadratic term aT^ due 
to the thermal excitation near the nodal region. Thus 
A7 can be written as: 
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FIG. 5: Field induced DOS at zero K of the overdoped sam- 
ple (x=0.22). The solid line is a theoretical curve A7 — 
1.995VH{1 - 0.0121//). The dotted line represents the best 
fit to the case at the unitary limit (refi— ). The inset shows 
the slope dA^/dT of the straight lines shown in Fig. 4 in low 
temperature region. The dashed line is a linear fit to the data 
at fields above one tesla. The intercept and the slope of the 
dashed line give rise to the pre- factors of the second and third 
terms in eq.7 



A7 = 7(iJ) - 7(0) = hVH + (62 - a)T + bcoreH (5) 

From eq.5 one can see that A7 depends on T through 
the second term, however the slope 62 — a is still field 
independent by definition. This clearly indicates that 
the Simon-Lee scaling law is still not enough to interpret 
the field dependent slope of A7 vs. T as shown in Fig. 4. 

Let us keep going, still based on eq.5, we propose that 
the core size effect may have a sizable influence on the 
total vortex quasiparticle excitations. This is actually 
reasonable since the vortex core with size ^ ex hvp/Ag 
grows up in the overdoped side due to smaller supercon- 
ducting gap value^^, where vp is the Fermi velocity and 
As is the superconducting gap. By taking the vortex 
core size (2^) into account, i.e., deducting the normal 
core area away from the Volovik term, one can rewrite 
A7 as 



A7 = {hVH + b2T) X (1 - i'/Rl) -aT + horcH (6) 

where ^ is the radius of the normal core, Ra is the 
outer radius of a single vortex where the supercurrent is 
flowing, thus R^ — (f)a/T^H. Reorganizing all terms in 
eq.6 leads to 



<^ 



ne 



A7 = 6iVi7x {l-^H) + {b2-a)T-b2^HT+b,oreH 

(7) 

One can see that the third term in eq.7 is just what we 

need for interpreting the difficulty as mentioned above, ft 



is necessary to recall that the core size correction is pro- 
portional to ^^, for example, it will be four times when 
^ doubles. Thus increase of ^ in the overdoped side will 
give sizable effect on the total DOS and core size effect 
should be considered. Next let us have a closer inspec- 
tion at the data and derive some parameters. At zero 
temperature, only the first term and the last term are 
left in eq.7. The values of A7(T = 0) are determined 
from the extrapolation of the linear lines in Fig. 4 to 
K and presented in Fig. 5. The data A'y{H,T = 0) is 
also determined by doing linear fit to the raw data C /T 
vs. T^ between 2K and AK, and then subtracted the 
zero temperature value 70. The results are quite close to 
each other by using these two different methods. The 
solid line in Fig. 5 is a fit to the data using the first 
term in eq.7 yielding 61 = 1.995 ± Q.M&mJK-'^T-^/'^ 
and TT^V^o = 0.012 ± 0.003 and thus S, = 28. 2A. The 
value ^ = 28A derived here is quite close to that found 
in Nernstffi and STM measurements^^ (20A for opti- 
mally doped Bi-2212 sample). We also tried to use the 
first term together with the last term to fit the data 
but find out that the contribution from the last term 
is extremely small. The first term here describes the 
zero temperature data very well, indicating the absence 
of a second component of order parameter such as id^y 
or is since otherwise the Fermi surface would be fully 
gapped and the Doppler shift had very weak effect on 
the quasiparticle excitations. The inset of Fig. 5 shows 
the field dependence of the slope of the linear part in 
Fig. 4. It is clear that the slope increases roughly lin- 
early with H above 1 T. This can be exactly anticipated 
by the second and third terms in eq.7. From the in- 
set of Fig. 5 one obtains a — 62 = 0.305mJmo/~^i^~'^ 
and 627r^^/0o = 0m238mJmol-^K-^T-^. By tak- 
ing ^ = 28.2A, we obtain the following values: a — 
0.501mJmol-^K-^ and 62 = 0.196mJmol-^K-^. The 
value of a = 0.501mJmol^^K^^ found here is quite close 
to the value obtained by fitting the zero field data to eq.2 
(0.465 mJ/molK^ see Table-I). Since the contribution 
from the core region (last term in eq.7) is negligible com- 
paring to the Volovik term, from eq.7 one understands 
that the failure of using the Simon-Lee scaling law in very 
overdoped sample is due to the core size effect. This is 
actually quite reasonable since the core size in the over- 
doped region grows up. Sofar we don't know yet whether 
the negligible contribution from the core region is because 
of the gapped feature within the core region as found in 
optimally and underdoped samples^iiSiSSiSLSSiSS, or it is 
naturally small comparing to the contributions from the 
Doppler shift effect of the surrounding superfluid. This 
casts a interesting issue for future STM measurement on 
the tunnelling spectrum within the vortex cores in very 
overdoped region. 

It is necessary to estimate how much of the field in- 
duced delocalized DOS is contributed by the impurity 
scattering in our present sample. At the unitary limit 
at zero energy, i.e., when T=0, Kiibert and Hirschfeld^ 
predict that the field induced relative DOS is S-f/^g = 
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FIG. 6: Scaling of the raw data [C{H) - C{0)]/TVH vs. 
T lyH for the overdoped sample based on the Simon-Lee scal- 
ing. Clearly no good scaling can be obtained. 



FIG. 7: Plot of the data C.^^i = {AC/T^ + a)/{l-n£,^H/(l)Q) 
vs. T/\/7l. It is clear that the data collapse onto one main 
branch which is expected by the theoretical expression with 
core size correction (eq.8). 



PiiH/P2)logiP2/H), where Pi = 0.322(Ao/r)i/2, Aq 
the gap maximum, T the impurity scattering rate, and 



P2 = 7r7Je2/2a2 



1. The dotted hne in Fig. 5 rep- 



resents the best fit of this relation to our data yielding 
r/Ao ~ 0.00039 (close to the clean limit). In addition, 
the value of Hc2 derived here is about 21.53T, which is 
too small for present sample. It is clear that the fit has a 
poor quality comparing to the better fit in the clean limit 
(solid line). Furthermore, the formulism considering the 
impurity effect does not predict a field dependent slope 
for the linear relation A7 vs. T in low temperature re- 
gion. Therefore together with the extremely small F/Aq 
found in present case, we believe that the field induced 
DOS in our sample comes mainly from the Dopplcr shift 
effect on supercurrent outside the cores. The residual 
linear term ( 70T ) of electronic specific heat will be dis- 
cussed separately in the following subsection. 

In order to show the inapplicability of Simon-Lee scal- 
ing law for the overdoped sample, we present the raw 
data [C{H) - C{0)]/tVH vs. T/y/H in Fig.6. Clearly 
the scaling looks very poor. From above discussion, we 
conclude that the failure of the Simon-Lee scaling law in 
the very overdoped region is due to the quite large vortex 
core size which needs to be corrected. 

The nice fit in Fig. 5 with only the first term of eq.7 
suggests that the core region has very small contribution 
to the DOS since otherwise the last term bcoreH should 
be sizeable. This implies that the low energy DOS inside 
the vortex core is very small. Based on this idea, we write 
a new scaling law as 



AC/T^ 



1 - Tren/cPo 






(8) 



One can use this equation to test the idea about the 
vortex core size correction. We thus present the data 
of (AC/r2 + a)/(l - 7r^2^/$o) vs. T/^/H in Fig.7 



with a = Q.501mJ/molK^ and tt^V'/'o = 0.0121 de- 
rived above. The data collapse on one branch and show 
good consistency with the expected theoretical curve. 
The slight scattering or deviation from the main scal- 
ing branch is due to the simple assumption made for the 
core size correction (1 — C^/i?^) and the rough estima- 
tion for a value. Worthy of noting is that to have this 
nice data collapsing and consistency with the theoretical 
curve we need to take bcore ~ 0, again showing a small 
contribution from the inner vortex core. The nice data 
collapsing using eq.8 suggests that the Simon-Lee scaling 
law can be reconciled by considering the vortex core size 
effect. It is interesting to note that the electronic ther- 
mal conductivity derived by Sun et ali^ is not consistent 
with the Volovik's expression in low temperature region 
for the overdoped sample, rather it shows a plateau when 
the field is high. However the H^^^ law is followed very 
well in low temperature region for the optimally doped 
sample. Our core-size correction picture may provide al- 
ternative interpretation to this discrepancy. 



C. Optimally doped sample (x=0.15) 

In order to have a comparison with the overdoped sam- 
ple, in this subsection we present the data from an op- 
timally doped one with x=0.15. The raw data of spe- 
cific heat for the optimally doped sample is shown in 
Fig. 8. The separation between each field can also be 
easily distinguished in low temperature region. Again 
here the curve at zero field extrapolates to a finite value 
(70) at K instead of zero. This will be discussed in 
the forthcoming subsection. It is found that the lin- 
ear behavior of A7 vs. T for the overdoped sample 
(shown in Fig. 4) is absent here. This may be due to 
the much smaller a value (see Table-I). We then check 
whether the d-wave scaling law is applicable here. If 




T' ( K' ) 



FIG. 8: Raw data of C{H)/T vs. T^ for the optimally doped 
sample. The inset shows the diamagnetic transition measured 
in the ZFC mode at 20Oe. 
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FIG. 11: Zero temperature specific heat Aj{H) of the opti- 
mally doped sample. The solid line is a theoretical expression 
A-y(H) = 0. 577^17 (mJ/molK^) which fits the data very 
well. The dashed line is a fit to the first term of eq.7 yielding 
a small and unreasonable negative value for tt^^ /cj)o. 
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FIG. 9: Scaling of the raw data [C{H) - C{Si)\lT\/li vs. 
T I^TH for the optimally doped sample (x = 0.15). The scaling 
looks rather good. 
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FIG. 10: Scaling of the raw data [C(H)-C(0)]/r2 vs. r/\/ff 
for the optimally doped sample based on the Simon-Lee scal- 
ing law. Clearly the scaling quality is quite good. 



the Volovik (Dopplcr shift) effect really dominates here, 
one can expect that C{H) - C(0) = T^y{T/^) - aT^, 



thus [C{H) - C(0)]/T2, or [C{H) - C{0)]/TVH should 



scale with T/^/H. In Fig. 9 we present the result of 
[C{H) - C{0)]/Ty/H vs. T/^/H. It is clear the scal- 
ing is rather good comparing to that of the overdoped 
sample [Fig. 6]. Here the value of A^/^/lI in the zero 
temperature limit gives the pre-factor A in the Volovik's 
relation C^i — AT\/H vi^hich is about 0.55 to 0.6 
mJmol~^ K~^T^^/^ . In Fig. 10 we present the Simon-Lee 
scaling in the way of [C{H) - C(0)]/T2 vs. T/^/H. One 
can see that the scaling is reasonably good. All data be- 
low about 10 K collapse onto one branch. We have been 
aware that Nohara et alJ^ successfully used the Simon- 
Lee scaling law to the overdoped sample x = 0.19, but 
failed for the optimally doped one. The failure of using 
Simon-Lee scaling law in Nohara's experiment for opti- 
mally doped sample is in contradiction with the reports 
from many other groupsiiiiSiil. This may be caused by 
the way that they used to derive A7. As stressed in pre- 
vious subsection, we use A7 = [Cjy ||c — Ch=o\/T instead 
of using A7 = [Cj:/||c — Chlc]/T to derive the field in- 
duced change of 7. The latter (as used by Nohara et al.) 
may inevitably involve the unknown DOS contributions 
from another kind of vortices (for example, Josephson 
vortices) when H J- C. For La — 214 system, since the 
Schottky anomaly is very week, it is not necessary to 
derive A7 in the second way. While Nohara et al.^^ ob- 
tained a relatively good scaling for the overdoped sample 
{x = 0.19). We would not comment on the validity of 
this successful scaling at a; = 0.19. One reason for the 
discrepancy between their results and our results may 
be from different doping levels: our sample {x = 0.22) 
is more overdoped and the vortex core size is certainly 
larger and needs to be corrected . 

The data Aj{H, T = 0) is determined by doing linear 
fit to the raw data C/T vs. T^ between 2K and iK, 



and then subtracted the zero temperature value 70. The 
results are shown in Fig.ll. We tried to fit the zero tem- 
perature data in Fig.ll to the first term in eq.7 in terms 
of core size correction ( shown by the solid line) , it turns 
out that the correction term 7r^^i?/0o is very small and 
negative, which is certainly unreasonable. This actually 
indicates that A7(T = 0) can be nicely fitted to the theo- 
retical expression A^{H) — 0. 577 ^/H{mJ/m,olK^). Us- 
ing C'yoi/T = k^n\/H/Hc2, we have kjn = 0.577^/5c2■ 
A similar value (0.49) was derived by Fisher et al>ii for 
ia-214 sample with X = 0.15. Taking i/c2 ~ lOOTiS^and 
k Ri 0.74ii, we have 7„ « 7.8mJ/molK^, which is very 
close to the reported values for optimally doped La — 21A 
sampleiii^. Worthy of noting here is that the field in- 
duced extra DOS at zero K can be nicely fitted with the 
Volovik's relation Cyoi « T^/H albeit the residual linear 
term 70 is quite large. This suggests that the residual 
linear term 70 observed commonly in curate supercon- 
ductors may originate from some other properties, such 
as inhomogeneity. ft may not be induced by the small 
scale impurity scattering since otherwise the H^'^ law 
should not be followed so well. Since both the Simon-Lee 
scaling law and the Volovik's ^fH are followed very well 
for optimally doped sample, the core size effect seems to 
be very week. 
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FIG. 12: AC susceptibility and DC magnetization of the 
underdoped sample with p=0.069. The bottom curve is the 
real part x' of the ac susceptibility, and the upper one is the 
imaginary part x" . By extrapolating the most steep transition 
portion of the real part of the AC susceptibility to the normal 
state background {x' = 0), the Tc = 12K is determined here. 
The inset shows the DC magnetization measured in the FC 
and ZFC processes. Below about 7.5 K, the magnetization 
measured in FC mode is rather stable. The slight temperature 
dependence of the magnetization measured in ZFC mode is 
induced by the flux penetration. 



D. Underdoped sample 

In this subsection the low temperature specific heat of 
underdoped La2-a;Sr^Cu04 {p = 0.069,0.075,0.09 and 
0.11) single crystals is reported in magnetic fields up to 12 
T. It is found that the Volovik's relation C^oi = ATH^/"^ 
is still satisfied in the zero temperature limit, but the pro- 
posed Simon-Lee scaling law, i.e., Cvoi/T"^ = f{T/\/H), 
is not followed so well, except for at very low tempera- 
tures (below about 3-4 K). 

Fig. 12 shows the temperature dependence of the AC 
susceptibility and DC magnetization of the underdoped 
sample La2~xSrxCuOi {p — 0.069). The transition tem- 
perature drops from about 12 K to 9 K after extract- 
ing some oxygen out of the sample (not shown here) 
by annealing the sample in Ar gas for 48 hrs. Then 
Tc keeps stable upon further annealing in Ar gas. The 
Tc is increased again when the sample is treated in flow- 
ing oxygen. The DC magnetization measured in the FC 
process shows a transition width of about 2.5K. Below 
about 7.bK the M(T) curve keeps flat. The magnetiza- 
tion measured in ZFC mode shows a slight increase with 
temperature induced by the easy flux penetration in the 
very underdoped region. Specific heat has been mea- 
sured in the FC mode as done for all other samples. This 
mode provides a vortex system which is close to equi- 
librium state and thus relatively uniform^. Presented 
in Fig. 13 is the specific heat C/T as a function of T^ 
at magnetic fields ranging from to 12T for this un- 
derdoped sample before annealing (estimated p=0.069). 
fn low temperature region the curves are rather linear 



showing that the major part is due to the phonon con- 
tribution Cph = I3T^ -f (5r^, and has no slight upturn 
in low temperature region due to the Schottky anomaly 
of free spins. The curve at zero field extrapolates to a 
finite value (70) at OK , again showing the existence of 
a residual linear term which will be discussed later. As 
mentioned before the phonon part is independent on the 
magnetic field, this allows to remove the phonon contri- 
bution by subtracting the C/T at a certain field with that 
at zero field, one has AC = C{H) - C(0) = C^oi - aT^ 
and AC/T^ = Cyoi/T'^ - a. The results after the sub- 
straction are shown in Fig. 14. One can see that the linear 
part with negative slope as appearing for the overdoped 
sample is absent here. This is understandable when a 
value (or aT^ term) is very small comparing to the field 
induced change of total specific heat. Therefore for this 
underdoped sample, no apparent T^ term at H = was 
observed which can be found easily in the overdoped 
LSCO sample. This is consistent with the data shown 
in Table-I and experimental results from other groups on 
LSCOiii4. The disappearance of this aT^ term was usu- 
ally interpreted as due to either the impurity scattering 
which smears out the nodal effect, or the small value of 
coefficient a of the T^ term. We will show that this is 
induced by much smaller 7„ value (thus smaller a) just 
above Tc- 

Next let us have a look at the field induced DOS at 
T=0 K. The data A-/{H,T = 0) is determined by doing 
linear fit to the raw data C/T vs. T^ between 2K and 
AK, and then subtracted the zero temperature value 70. 
The results are shown in Fig. 15. In order to compare with 
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FIG. 13: Raw data of C{H)/T vs. T^ for the underdoped 
sample p=0.069. One can see that the field induced change of 
specific heat becomes much smaller than that of the optimally 
and overdoped sample. 



FIG. 15: Field induced change of 7 at zero K for the un- 
derdoped sample with p=0.069. The solid line and dashed 
line are fits to the theoretical relation of Volovik effect at the 
clean limit and the impurity scattering at the unitary limit 
respectively. 
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FIG. 14: The subtracted data [C(H)-C(0)]/T vs. T for the 
underdoped sample [p = 0.069). It is clear that no linear 
part with negative slope of A7 vs. T as appearing for the 
overdoped sample can be observed here. This may be induced 
by the much smaller a value. 



FIG. 16: Scaling of the raw data I^CIT^JH vs. T/V-ff for 
the sample p = 0.069 based on the Simon-Lee scaling. Clearly 
no good scaling can be obtained, except for at very low tem- 
peratures. 



the theoretical predictions, the increase in ^{H,T = 0) 
was fitted with Aj{H,T = 0) = AH^ , and the value 
of B is 0.52 and A is about 0.28. The value B derived 
here from free fitting is very close to 0.5 as predicted by 
the Volovik theorjii^ which may manifest the existence 
of hne node in the gap function. We can also fix _B = 0.5 
and find out that A = 0.282mJ/molK^T^/^. This is also 
compatible with the results of other groups^^. For the 
zero temperature data we also considered the core size 
correction, i.e., tried to use the first term of eq.7 to fit the 
zero temperature data. But it gives rise to a small and 
negative value of tt^,'^ /(f>o which is certainly unreasonable. 
For the underdoped sample, we used the Simon-Lee 
scaling law to scale our data. The results of [C{H) — 
C{0)]/Ty/H vs.T/Vh are shown in Fig.l6. The data 



fan out showing a poor scaling quality. Clearly the data 
cannot be scaled using the Simon-Lee scaling law except 
for at very low temperatures. We plot also the data of 
[C{H) - C(0)]/r2 = Cdos/T^ - a vs. T/y/H in Fig.l7, 
one can again see the poor scaling in wide temperature 
region. The Simon-Lee scaling has been applied to all 
samples investigated in this work (p= 0.069, 0.075, 0.09, 
0.11, 0.15, 0.22). It is easy to find that the scaling qual- 
ity becomes better and better when the doping concen- 
tration is increased from 0.069 to 0.15. One can even 
see the gradual change among these underdoped sam- 
ples (p=0.069, 0.075, 0.09, 0.11): the scaling curves fan 
out hke that in Fig. 17 for samples with p = 0.069, 0.075, 
but the scaling pattern becomes narrower towards higher 
doping. The scaling behavior are shown in Figs. 18-20 
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FIG. 17: Scaling of the raw data AC/r^ = [C{H)-C{0)]/T'^ 
vs. T/VH for the sample p = 0.069 based on the Simon-Lee 
scaling. Good scaling can be found only at very small values 
of T/VH. 



FIG. 19: Scaling of the raw data AC/T^ = [CiH)-C{0)]/T'^ 
vs. T/VH for the sample p = 0.09 based on the Simon-Lee 
scaling. 
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FIG. 18: Scaling of the raw data AC/r^ = [C(H) -C(0)]/r2 
vs. T/\/Tl for the sample p — 0.075 based on the Simon-Lee 
scaling. Good scaling can be found only at very small values 
of T/VH- 



FIG. 20: Scaling of the raw data AC/r^ = [C{H)-CiO)]/T^ 
vs. T/y/H for the sample p = 0.11 based on the Simon-Lee 
scaling. Now the fanning-out of the scaling curves are strongly 
constrained showing a better scaling behavior. 



for samples with p — 0.075, p — 0.09 and p = 0.11. A 
clear trend for a better scaling at a higher doping can be 
easily seen here. 

There are several possibilities for the failure of using 
Simon-Lee scaling law in underdoped region. One possi- 
bility is due to the impurity scattering effect as suggested 
by Kitbert et al.~ . Thus we use the dirty limit formula 
^{H) = 7(0)[1 + D{H/Hc2)ln{Hc2/H)] to fit the data at 
zero K, where D w A/32r. For the simplicity we show 
here only the fit to the data of the sample p ~ 0.069. It is 
found that the data can also be roughly fitted by the re- 
lation with impurity scattering (as shown by the dashed 
line in Fig. 15). The obtained results for the sample with 
p = 0.069 are: i/c2 = 45.6T, 7(0) = 4.03mJ/mol K^, 
r/A = 0.046. Thus it seems that one cannot rule out 
the possibility of impurity scattering to play a dominant 



role in the field induced change of 7 in very underdoped 
sample. However this speculation cannot interpret the 
nice \/H dependence of the field induced DOS at zero K 
as shown in Fig. 15. Worthy of noting is that the dirty 
limit formula of Kiibert et al.^° is more flexible to fit to 
the data than the simple VH relation. One needs to seek 
an alternative way to clarify this discrepancy. 

The second possibility is the core size effect as appear- 
ing in the overdoped sample. We then try to scale the 
data by using eq.8 and leaving both a and ir^"^ /(J)q as 
free fitting parameters. Unfortunately no good scaling 
can be found by choosing any values for a and 7r^^/0o- 
This is in consistent with the fact that an unreasonable 
negative value for tt£,'^ /(f>o is obtained if we fit the zero 
temperature data in Fig. 15 to the first term of eq.7. Both 
indicate that the failure of Simon-Lee scaling law here is 
not due to the core size effect. One may argue that the 
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data is scalable with only a very narrow scaling region of 
T l\fli, for example, from Fig. 16-18, the scalable region 
is about T/\/F < l.SKT-o^ This is of course possible 
since we don't know the precise value for many parame- 
ters. However we can have a rough estimation to check 
whether this is rea sonable . Provided the scalable region 
is T|^Jli < Tc/y^Hc2iO) = 1.5, inputting T^ = UK, 
one has i?c2(0)=64 T which seems too big for this very 
underdoped sample. 

Another possible reason for the failure of the scaling 
law is that the sample is in the underdoped region with a 
pseudogap in the normal state. When the sample is in the 
mixed state, some competing or coexisting order such as 
short range antiferromagnetic orderi^iSLSLSSi^SiiMi, or 
the SDW orden^, or a d-density wave (DDW) ordei^ 
is enhanced, and this newly generated or enhanced or- 
der will certainly give contribution to the total specific 
heat. For example, for 2D AF correlation, it is known 
that Caf oc T^. Therefore qualitatively the failure of 
the Simon-Lee scaling law in underdoped region can be 
understood in the following way. By increasing the mag- 
netic field, a second order is generated or enhanced within 
the vortex core and nearby regions (about lOOA). On 
one hand this region is gapped leading to the decrease 
of the total DOS at fermi level simply by reducing the 
region where the supercurrent can flow. On the other 
hand the newly generated AF or SDW or DDW region 
will contribute a new term to the total specific heat due 
to spin or other type excitation. The relevant competing 
order under a magnetic field, according to both neutron 
scatteringi2iS2i and NMR measurement^, may be the AF 
order. STM measurement by Hoffman et al.^'' indicates 
a checkerboard like modulation with periodicity of 4a of 
the LDOS. This was regarded as the direct observation 
of the strong electronic correlation with the underlying 
competing order which was predicted by many theoreti- 
cal work'^Si39i40,4ii45,46,47,48^ rpj^j^ qualitative picture calls 

for further detailed analysis and evidence from other ex- 
periments. Since the heat capacity from the newly gener- 
ated or enhanced second order has a temperature depen- 
dence of T"^ with e > 1, at zero temperature the specific 
heat from this term is zero, thus the H^'^ law from the 
Doppler shift of the d-wave superconductivity is restored. 
This may be the reason for that the zero-temperature 
data follows the ff^/^ law but the data at finite temper- 
atures do not satisfy the Simon-Lee scaling very well. 



E. The residual linear term 70 

Almost in all cuprate superconductors, a residual lin- 
ear term of electronic specific heat 70 has been observed 
in the low temperature limit T —^0, even in the best sam- 
ples up to date. In YBa2Cu2,0'j single crystals, Moler et 
dlr- found that 70 is larger for the twinned samples than 
the detwinned ones. Meanwhile they further found that 
7o increases when the sample becomes more underdoped. 
Clearly one can conclude that 70 is related to the impu- 
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FIG. 21; Field dependence of the field induced extra 7 nor- 
malized by the value at 12 T at zero K. The solid line rep- 
resents the theoretical curve H^'^. It is found that the data 
from different samples at different doping levels follow the 
H^'^ law reasonably well 
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FIG. 22: The doping dependence of 70 of La2-xSrxCu04, 

crystals. It is clear that the minimum 70 value is found in 

the region around p =0.11 to 0.125 . 

rities or disorders in the samples. While quite surpris- 
ingly, for many samples with quite different 70 values, it 
is found that the zero temperature data can be expressed 
as ^{H,T = 0) = 7o -h A\/H, showing evidence for d- 
wave pairing. This may suggest that 70 is not mainly 
induced by the impurity scattering, since otherwise the 
field induced extra DOS should not follow the relation 
AC{H,T = 0)/T = AVh so weU. In Fig.21, we present 
the field dependence of A7 = [C{H,T = 0) - C(0,T = 
0)]/T normalized to the value for each sample at about 
12 T. Meanwhile we show the H^^^ law by the solid line. 
One can see that for almost all samples, the field induced 
extra DOS at zero K follows the H^'^ relation reason- 
ably well despite the 70 value highly disperses. This fea- 
ture was also discovered by Chen et al>i^ on three typical 
samples (x = 0.10,0.16 and 0.22). Nohara et alii^ mea- 
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sured three single crystals (x — 0.10, 0.16 and 0.19) and 
empirically found that the optimally doped sample (x 
= 0.16) has the lowest value of 70 (70=2.8, 1.5 and 2.2 
mJmol~'^K~'^ for x=0.10, 0.16 and 0.19 respectively ob- 
tained from Fig.l of Ref.^"^). Chen et al.^"' found the simi- 
lar behavior among three samples with x = 0.10, 0.16 and 
0.22 (70=1.49, 0.7 and 1.41 mJmol~^K-^ respectively). 
This raises the question of the origin of this residual linear 
term and its correlation with the field induced quasipar- 
ticle DOS. As mentioned before, if the field induced DOS 
is related to the impurity scattering, another relation^S 
(57/70 = Pi{H / P2)log{P2 / H) is expected. This is some 
time contradicting to the experimental result (see dotted 
line in the main panel of Fig. 5). In addition, the yH 
dependence of the field induced change of 7 is certainly 
not obtained by accident, since it is found on different 
samples from different groups, even on poly-crystalline 
samples^*. From the point of view of chemistry it is 
not true that the optimally doped sample is the clean- 
est one since in most cases the underdoped samples can 
be more easily grown with high quality. In this sense 
the residual linear term may be related to some other 
properties rather than the impurity scattering. In Fig. 22 
the 7o values for different single crystals measured in our 
experiment from underdoped to overdoped are shown. 
The value of 70 is obtained by fitting zero field data 
to eq.2 (see Table-I). It is clear that the minimum 70 
is found in the region around 0.11 or 0.125. The value 
for 7o found from our data are more close to the data 
of Nohara et al,^'^ on single crystals, but clearly higher 
than that obtained on polycrystalline smplesi^. So far 
we don't know the reason for this discrepancy. For our 
extremely underdoped sample (x = 0.063) investigated 
here, although the data at finite temperatures cannot be 
treated with the Simon-Lee scaling law, the data in the 
low temperature limit T — > can be however still nicely 
expressed by 7(i:f, T = 0) = 70 + AvH, even the ab- 
solute increase of 'y{H, T = 0) by field is much smaller 
than 7o. Therefore it is reasonable to conclude that the 
field induced part is mainly contributed by the Doppler 
shift effect on the supercurrent outside the vortex cores, 
while the residual linear term 70 is mainly contributed 
by some small normal regions which weakly depends on 
the magnetic field. Similar explanations to the origin 
of 7o have been suggested for many times in the past^. 
This may be understood in the following way. In un- 
derdoped Bi — 2212 single crystals, scanning-tunneling- 
microscopic (STM) measurement indeed reveal a mixture 
of superconducting regions with sharp quasiparticle co- 
herent peaks on the tunnelling spectrum, and the non- 
superconducting regions with pseudogap-like tunneling 
spectrum^''. In the overdoped side, the tiny normal cores 
as proposed in the Swiss cheese model^^, or the mesoss- 
copic normal regions suggested by Fukuzumi et al.^— and 
Wen et al.^"^ will contribute a residual term 70 which does 
not show an apparent increase with the field. As pro- 
posed by Fukuzumi et aliS that the dome-like electronic 
phase diagram may be formed by the mixture of three 



phases: anti-ferromagnetic phase in the extremely under- 
doped region, a d-wave superconducting region with the 
robust superconductivity near the optimal doping point 
and a non-superconducting Fermi liquid in the overdoped 
region. According to this simple picture the 70 should in- 
crease in the underdoped and overdoped region, which is 
just the case as shown by the data in Fig. 22. Therefore we 
would argue that the residual linear term may be mainly 
contributed by some non-superconducting regions due to 
phase separation, either chemical or electronic in origin. 
This interesting argument needs certainly to be further 
checked with data obtained by different techniques on 
different systems. 



F. Field induced reduction of specific heat in high 
temperature region 

In above analysis, we concentrate on the data below 10 
K ( below 6 K for the very underdoped sample). This is 
also the temperature region that most of the low temper- 
ature specific heat data was reported in the literatures. 
Now we report another phenomenon: field induced re- 
duction of specific heat in the mixed state. In Fig. 23 we 
present the temperature dependence of the field (12T) 
induced change of 7 for three typical samples analyzed 
above, here A7 = [C(12T) - C(OT)]/T. Although the 
data are strongly scattered one can still see that: (1) 
The field induced change A7 becomes negative at about 
0.5 — 0.7Tc; (2) The curves have a similar shape: A7 is 
positive in low temperature region, then it becomes neg- 
ative and finally comes back to zero in high temperature 
region (near Tc for optimal and overdoped sample). For 
the overdoped sample, the A7 keeps negative above 0.5Tc 
until Tc at which A7 suddenly goes back to zero. For the 
optimally doped one, the A7 is negative above about 
0.7Tc up to the highest temperature we measured here 
(30 K). However for the underdoped sample, it shows 
that the A7 keeps negative until 1.5Tc. Similar data were 
obtained by Fisher et al.^^ on samples with x = 0.15. 
Our data near Tc is more scattered since our setup can 
only measure samples with maximum mass of 50 mg. 
This feature, namely the negative A7 in high temper- 
ature region is a consequence of entropy consideration, 
which has been observed in all types of superconductors. 
In low temperature region, when a magnetic field is ap- 
plied, vortices will be generated leading to higher DOS 
near Fermi surface, so that A7 = ^{H) — 7(0) is posi- 
tive. When the temperature is increased, to satisfy the 
field independent entropy above Tc, in a certain region 
below Tc, A7 should be negative. The most interest- 
ing point is for the underdoped sample here, even above 
Tc one clearly sees a magnetic field induced change of 
entropy. This implies an abnormal normal state which 
is far from a conventional metal. For a conventional s- 
wave superconductor, the field induced change of 7 can 
be negative near Tc- It is difficult to understand the field 
induced reduction of specific heat well below Tc, since the 
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normal core region always gives rise to a higher DOS of 
quasiparticle. Outside the vortex core the DOS is almost 
negligible. However this field induced reduction of spe- 
cific heat well below Tc is found to be a general feature 
of all LSCO crystals we investigated so far. This may 
be related to the intrinsic properties of cuprate super- 
conductors. In a d-wave superconductor, theoretically it 
is predicted that there is a ZBCP within the vortex core 
which should also contribute a quite high DOSSSiS^ in the 
mixed state. Besides, a high DOS will be generated by 
the Doppler shift effect of the supercurrent surrounding 
the vortex core. Normally the sum of these two terms 
are larger than the zero-field term Cdos = a^^, leading 
to a field induced enhancement of DOS in low tempera- 
ture region. When the temperature is high, the Doppler 
shift effect will be smeared out by the strong thermal ex- 
citation and finally A7 becomes negative. As far as we 
know, no quantitative theoretical expression about A7 
has been reported so far for a d-wave superconductor in 
wide temperature region. We cannot have a quantitative 
understanding to our data. However, this field induced 
reduction of specific heat well below Tc may be under- 
stood as due to the anomalous feature of vortex core 
state, i.e., a gapped vortex core as seen by the STM1S&, 
or based on the assumption that the contributions from 
the core region is much smaller than the outside region 
where either the Doppler shift or the strong thermal ex- 
citation dominates. Actually the Simon-Lee scaling law 
becomes a C^oi oc T^ relation in high temperature re- 
gion. In this case the quasiparticle excitation outside the 
vortex core is almost the same {aT'^) with or without 
applying a magnetic field. However since the vortex core 
region is gapped or contributes negligible part to the to- 
tal DOS, one needs to take the core region away from 
the total area in calculating A7, naturally leading to a 
negative value of A7. 



IV. DISCUSSION 

In low temperature region, our analysis indicates that 
the field induced quasiparticle DOS can be well described 
by Volovik's theory or Simon-Lee scaling law although a 
correction due to the core size effect is needed for the 
overdoped sample. This means that the prerequisite for 
the theory, i.e., the dx2-y2 pairing symmetry is well sat- 
isfied. Therefore it naturally rules out the presence of 
a second order parameter like id^y or is either due to 
overdoping^ or due to the field effect-'* in all samples in- 
vestigated here. Meanwhile, for the overdoped sample, 
another interesting phenomenon is that the vortex core 
region contributes very little (at least much smaller than 
that induced by the Doppler shift if the super-current 
would flow in the same area) to the total DOS. We have 
also tried to analyze the data of the optimally doped 
and underdoped sample in the way as that for overdoped 
one, for example to fit the data in Fig. 11 and Fig. 15 to 
the first term in eq.7. It turns out however that the 
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FIG. 23: Temperature dependence of the field (12T) induced 
change of A7 = [C(12T) - C{OT)]/T vs. T for three typical 
samples (x=0.22, overdoped: filled circles; x=0.15, optimally 
doped: diamonds; x=0.069, underdoped: triangles ). The 
horizontal axis is normalized to Tc of each sample. In high 
temperature regime below Tc, the field induced change of DOS 
becomes negative. 



correction term tt^'^H/^q derived is small and negative 
which is unreasonable. For the optimally doped sample, 
it is quite easy to understand since the vortex core be- 
comes very small. However for the underdoped sample, 
it is quite hard to understand since the core size tends 
to grow up too2^. The negligible contribution from the 
vortex core region may suggest that the ZBCP is ab- 
sent within the cores, even in the overdoped region. This 
suggestion inferred from the specific heat measurement 
about the ZBCP within the vortex core is consistent with 
the tunnelling results2^iSS42SiSL2& and certainly clears up 
the concerns about the surface conditions in the STM 
measurement. Recent results from NMR also show the 
absence of a ZBCP inside the vortex core^i. In this sense 
our data together with the earlier NMR data present a 
bulk evidence for an anomalous vortex core. Interest- 
ingly it is widely perceived that the normal state in over- 
doped region shows a Fermi liquid behavior even when 
the superconductivity is completely suppressed^. If this 
is the case the mean-field frame of BdG theory based 
on the conventional d-wave superconductivity seems not 
enough to interpret the anomalous vortex core state in 
HTS. For the underdoped sample, the Simon-Lee scahng 
fails except for in very low temperature region. This is 
interpreted as due to the presence of a second (gapped) 
order like AF or SDW or DDW within and nearby the 
vortex core. However one needs more theoretical and ex- 
perimental efforts to show the justice of this argument. 

By fitting the field induced extra DOS at zero tem- 
perature to the relation A7 = AH^^^, we obtained 
the pre-factor A in wide doping regime, where A = 
0.747„/V^c2^- The results are presented in Fig. 24. It is 
seen that the A-value increases with the doping concen- 
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FIG. 24: The doping dependence of the pre-factor A in 7(T = 
0) = 70 + AH^,B « 0.5. It is evident that the A-value 
increases with the hole concentration monotonously. The data 
measured on poly crystalline samples are somewhat smaller, 
which is perhaps induced by the random orientation of the 
grains. For some grains the field is not parallel to the c — axis 
leading to a smaller contribution to the field induced change 
of 7. 



tration monotonously. This can be understood in the fol- 
lowing way: by increasing doping the normal state value 
7„ will increase^, the Hc2 will drop down (at least it 
is the case in the overdoped region). Therefore A-value 
will increase monotonously in the overdoped side. One 
can see from the data that the A-value keeps almost con- 
stant in the extremely underdoped region, which means 
that 7„ and Hc2 should both decrease with underdop- 
ing. This indicates that the Hc2 becomes smaller and 
the coherence length S, becomes larger towards more un- 
derdoping. This is consistent with the recent conclusion 
drawn by Wen et ali^ by analyzing the data about the 
low temperature flux dynamics. This conclusion about 
the coherence length calls for a direct check to the vortex 
core size by using scanning-tunncUing-microscopy in the 
future. 



V. CONCLUDING REMARKS 

In conclusion, the field induced change of the electronic 
specific in mixed state of a series La2-xSrxCuOi sin- 



gle crystals has been measured and extensively analyzed. 
It is found that the field induced DOS of the optimally 
doped sample fits the predicted Simon-Lee scaling law 
for a d-wave superconductor very well, while deviations 
have been found for the overdoped sample. However, it 
is reconciled for the overdoped sample if one considers 
the core size effect provided the contribution from the 
inner vortex core is small comparing to that due to the 
Doppler shift in the same area. The Simon-Lee scaling 
law is applicable in the underdoped region only in very 
low temperature region. We attribute this to the ap- 
pearance of a second competing order (like AF or SDW 
or DDW) within and nearby the vortex core. The negli- 
gible contribution from the vortex core region may sug- 
gest the absence of the ZBCP in the vortex core, even in 
the overdoped region, although it is expected by the Bo- 
goliubov de-Gennes theory for a d-wave superconductor. 
Finally we present the doping dependence of the residual 
linear term 70 commonly observed in cupratc supercon- 
ductors. It is argued that this linear term may be related 
to inhomogeneity (either electronic or chemical), rather 
than be simply explained as due to the small scale im- 
purity scattering as usually thought. This conclusion is 
made because the field induced extra DOS at zero tem- 
perature follows the Volovik's \fH law reasonably well 
in all doping regime. It is hard to believe that this nice 
consistency is obtained by accident. Our results gen- 
erally conclude a d-wave pairing symmetry for the hole 
doped La2~xSrxCu04^ samples, although some compet- 
ing orders may co-exist with the superconductivity, and 
an anomalous feature (missing of the ZBCP) may appear 
within the vortex core. 
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